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4.1  INTRODUCTION 


The  performance  of  an  aircraft  can  adequately  be  described  by  assuming 
the  aircraft  is  a  point  mass  concentrated  at  the  aircraft's  center  of  gravity 
(eg).  The  flying  qualities  of  an  aircraft,  on  the  other  hand,  cannot  be 
described  in  such  a  simple  manner.  The  flying  qualities  of  an  aircraft  must, 
instead,  be  described  analytically  as  motions  of  the  aircraft's  eg  as  well  as 
motions  of  the  airframe  about  the  eg,  both  of  which  are  caused  by  aerodynamic, 
thrust  and  other  forces  and  moments.  In  addition,  the  aircraft  must  be 
considered  a  three  dimensional  body  and  not  a  point  mass. 

The  applied  forces  and  moments  on  the  aircraft  and  the  resulting  response 
of  the  aircraft  are  traditionally  described  by  a  set  of  equations  known  as  the 
aircraft  equations  of  motion.  This  chapter  presents  the  form  of  the  aircraft 
equations  of  motion  used  in  the  Flying  Qualities  phase  of  the  USAF  Test  Pilot 
School  curriculum. 

The  theory  presented  in  this  chapter  incorporates  certain  simplifying 
assumptions  to  make  the  main  elements  of  the  subject  clearer.  The  equations 
that  will  be  developed  are  not  as  rigorous  (and  complicated)  as  those  used  for 
design  of  modern  aircraft,  but  the  basic  method  is  valid  and  will  provide 
analysis  techniques  that  are  accurate  enough  to  gain  an  insight  into  aircraft 
flying  qualities.  With  the  aid  of  high  speed  computers  the  aircraft 
designers'  more  rigorous  theoretical  calculations,  modified  by  data  obtained 
from  the  wind  tunnel,  can  often  give  results  which  closely  predict  aircraft 
flying  qualities.  This  is  of  substantial  benefit  in  the  development  cycle  of 
new  aircraft. 

4.2  OVERVIEW 

An  aircraft  has  six  degrees  of  freedom  (if  it  is  assumed  to  be  rigid), 
which  means  is  has  six  paths  it  is  free  to  follow:  it  can  move  forward, 
sideways,  and  down;  and  it  can  rotate  about  its  axes  with  yaw,  pitch,  and 
roll.  In  order  to  describe  the  state  of  a  system  that  has  six  degrees  of 
freedom,  values  for  six  variables  (unknowns)  are  necessary.  To  solve  for 
these  six  unknowns,  six  simultaneous  equations  are  necessary.  For  an 
aircraft,  these  are  known  as  the  aircraft  equations  of  motion. 

The  full  aircraft  equations  of  motion  (given  in  sections  4.10  and  4.11) 
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reflect  a  rather  complicated  relationship  between  the  forces  and  moments  on 
the  aircraft,  and  the  resulting  aircraft  motion.  The  derivation  of  the 
equations,  however,  follows  a  very  simple  pattern  starting  from  Newton's 
second  law  for  translational  and  rotational  motions. 

Newton's  second  law  for  translational  motions  is 
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where  F  is  the  sum  of  the  externally  applied  forces  and  mV  is  linear  momentum. 
Newton's  second  law  for  rotational  motions  is 
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where  G  is  the  sum  of  the  externally  applied  moments  and  H  is  angular  momentum. 

F  and  G  are  both  vector  quantities  which  can  each  be  represented  by  three 
component  equations  (corresponding  to  three  dimensional  space).  The 
translational  equation,  therefore,  describes  the  aircraft  with  respect  to  its 
three  translational  degrees  of  freedom,  while  the  rotational  equation 
describes  the  aircraft  with  respect  to  its  three  rotational  degrees  of 
freedom.  Newton's  second  law,  therefore,  yields  six  equations  for  the  six 
degrees  of  freedom  of  a  rigid  body. 
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The  Left-Hand  Side  (LHS)  of  these  equations  represent  the  applied  forces 
and  moments  on  the  aircraft  while  the  Right-Hand  Side  (RHS)  stands  for  the 
aircraft's  response  to  these  forces  and  moments.  Small  perturbation  theory 
will  be  used  to  linearize  these  equations  so  they  can  be  solved.  This  will 
also  yield  terms  known  as  stability  derivatives  which  indicate  the  influence 
of  various  aircraft  characteristics  on  the  resulting  aircraft  motions,  and  are 
useful  in  coirparing  aircraft,  calculating  MIL-SPEC  requirements,  etc.  The 
equations  will  also  be  used  to  derive  aircraft  transfer  functions  which  will 
be  a  fundamental  part  of  the  mathematical  modeling  of  the  aircraft  and  its 
control  system  in  later  chapters. 

Prior  to  beginning  the  derivation  of  the  aircraft  equations  of  motion,  a 
discussion  will  be  presented  of  aircraft  sign  conventions  and  coordinate 
systems.  There  will  also  be  many  abbreviations  and  symbols  used  during  the 
derivation  of  the  equations  of  motion  and  in  subsequent  flying  qualities 
chapters.  A  list  of  abbreviations  and  symbols  is  given  in  the  last  subsection 
of  this  chapter. 

4.3  SIGN  CONVENTIONS 


The  sign  convention  used  at  the  Test  Pilot  School  and  the  Flight  Test 
Center  defines  a  positive  control  movement  or  deflection  as  one  that  causes  a 
positive  aircraft  movement  (right  yaw,  pitch  up,  right  roll).  Figure  4.1 
shows  this  sign  convention  and  gives  the  positive  directions  for  many  of  the 
variables  that  appear  in  the  equations  of  motion.  The  TPS  (AFFTC)  sign 
convention  does  not  conform  to  the  convention  used  by  NASA  and  some  reference 
textbooks.  Sign  conventions  will  be  discussed  in  greater  detail  in  the  flight 
control  systems  chapter. 

4.4  COORDINATE  SYSTEMS 

There  are  many  coordinate  systems  that  are  useful  in  the  analysis  of 
vehicle  motion.  We  will  be  concerned  with  three  of  these  coordinate  systems: 
inertial,  earth  axis,  and  vehicle  axis.  According  to  convention,  all 
coordinate  systems  used  will  be  right-hand  orthogonal. 
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FIGURE  4.1.  VEHICLE  FIXED  AXIS  SYSTEM 
AND  NOTATION 


EARTH 
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4.4.1  Inertial  Coordinate  System 

An  inertial  coordinate  system  is  defined  as  a  system  in  which  Newton's 
second  law  is  valid.  The  equations  of  motion  must,  therefore,  be  determined 
in  an  inertial  coordinate  system.  Another  way  of  defining  the  inertial 
coordinate  system  is  to  assume  it  is  an  axis  system  fixed  in  space  that  has  no 
relative  motion  (Figure  4.2). 

Experience  with  physical  observations  can  be  used  to  determine  whether  a 
particular  reference  system  can  properly  be  assumed  to  be  an  inertial 
coordinate  system  for  the  application  of  Newton's  laws  to  a  particular 
problem.  For  space  dynamics  in  our  solar  system,  the  sun  axis  system  is  a 
sufficient  approximation  for  an  inertial  system.  For  aircraft,  the  earth  axis 
system  is  usually  a  sufficient  approximation  for  an  inertial  coordinate  system. 

4.4.2  Earth  Axis  System 

There  are  two  earth  axis  systems,  the  fixed  and  the  moving.  Both  will  be 
referred  to  with  the  letters  XYZ  for  the  three  coordinate  axes.  An  example  of 
a  moving  earth  axis  system  is  an  inertial  navigation  platform.  An  example  of 
a  fixed  earth  axis  system  is  a  radar  site  (Figure  4.3). 
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FIGURE  4.3.  THE  EARTH  AXIS  SYSTEMS 
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In  both  earth  axis  systems,  the  Z-axis  points  toward  the  center  of  the 
earth  along  the  gravitational  vector,  g.  The  XY-plane  is  parallel  to  the 
local  horizontal  while  the  orientation  of  the  X-axis  is  arbitrarily  defined 
(often  defined  as  North).  The  two  earth  axis  systems  are  distinguished  by  the 
location  of  their  origins.  The  origin  of  the  fixed  system  is  usually  taken  as 
an  arbitrary  location  on  the  earth's  surface.  The  origin  of  the  moving  system 
is  usually  taken  as  the  vehicle's  eg.  What  distinguishes  the  moving  earth 
axis  system  from  the  vehicle  axis  system  discussed  in  the  next  subsection  is 
that  the  moving  earth  axes  are  not  fixed  in  orientation  with  respect  to  the 
vehicle.  They  are  instead  fixed  with  respect  to  local  vertical.  In  the  rest 
of  this  chapter,  the  XYZ  (upper  case)  system  will  be  assumed  to  be  the  fixed 
earth  axis  system  unless  otherwise  noted. 

4.4.3  Vehicle  Axis  Systems 

These  coordinate  systems  have  origins  fixed  to  the  vehicle,  and  axes 
defined  with  respect  to  the  vehicle.  There  are  many  different  types,  four  of 
which  are  commonly  used  for  describing  aircraft  motion:  the  body  axis  system, 
the  stability  axis  system,  the  principal  axis  system,  and  the  wind  axis 
system.  The  body  and  the  stability  axis  systems  are  the  only  two  that  will  be 
used  during  this  course. 

4. 4. 3.1  Body  Axis  System.  The  body  axis  system  (Figure  4.4)  is  the  most 
general  kind  of  axis  system  in  which  the  origin  and  axes  are  fixed  to  a  rigid 
body.  The  use  of  axes  fixed  to  the  vehicle  ensures  that  the  moments  and 
products  of  inertia  in  the  equations  of  motion  are  constant,  to  the  extent 
that  mass  can  also  be  considered  constant,  and  that  aerodynamic  forces  and 
moments  depend  only  upon  the  relative  velocity  orientation  angles  a  and  0. 
The  body  fixed  axis  system  is  also  the  natural  frame  of  reference  for  most 
vehicle-borne  observations  and  measurements  of  the  vehicle's  motion  and  will 
be  referred  to  as  the  xyz  (lower  case)  system. 

In  the  body  axis  system  the  unit  vectors  are  I,  j  and  k  with  origins  at 
the  vehicle  eg.  The  positive  x-axis  points  forward  along  a  vehicle  horizontal 
reference  line  with  the  positive  y-axis  out  the  right  wing.  The  positive 
z-axis  points  downward  out  the  bottom  of  the  vehicle,  usually  such  that  the 
xz-plane  is  the  vehicle  plane  of  symmetry. 
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FIGURE  4.4.  BODY  AXIS  SYSTEM 

4. 4. 3. 2  Stability  Axis  System.  Stability  axes  are  specialized  body  axes 
(Figure  4.5)  in  which  the  orientation  of  the  vehicle  axis  system  is  determined 
by  the  equilibrium  flight  condition.  The  xg-axis  is  selected  to  be  coincident 
with  the  relative  wind  at  the  start  of  the  motion.  This  initial  alignment 
does  not  alter  the  body-fixed  nature  of  the  axis  system;  however,  the 
alignment  of  the  axis  system  with  respect  to  the  body  changes  as  a  function  of 
the  equilibrium  condition.  If  the  reference  flight  condition  is  not 


FIGURE  4.5.  STABILITY  AXIS  SYSTEM 
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symmetric,  i.e.  with  sideslip,  then  the  xs-axis  is  chosen  to  lie  on  the 
projection  of  the  true  velocity  (VT )  in  the  plane  of  symmetry,  with  zs 
also  in  the  plane  of  symmetry.  The  moment  of  inertia  and  product  of  inertia 
terms  vary  for  each  equilibrium  flight  condition.  They  are  assumed  constant, 
however,  in  the  equations  of  motion. 

4. 4. 3. 3  Principal  Axes.  These  are  a  special  set  of  body  axes  aligned  with 
the  principal  axes  of  the  vehicle  and  are  used  for  certain  applications. 
Principal  axes  are  defined  as  those  axes  where  all  of  the  products  of  inertia 
are  reduced  to  zero.  The  equations  of  motion  are  thus  greatly  simplified,  but 
it  is  difficult  to  accurately  describe  the  aircraft  motion  in  this  system. 

4. 4. 3. 4  Wind  Axes.  The  wind  axes  use  the  vehicle  translational  velocity  as 
the  reference  for  the  axis  system.  Wind  axes  are  thus  oriented  with  respect 
to  the  flight  path  of  the  vehicle,  i.e.,  with  respect  to  the  relative  wind, 
V  If  the  reference  flight  condition  is  symmetric,  i.e.,  VT  lies  in  the 
vehicle  plane  of  symmetry,  then  the  wind  axes  coincide  with  the  stability 
axes.  The  wind  axes  depart  from  the  stability  axes,  moving  with  the  relative 
wind,  when  sideslip  is  present. 

The  relationship  between  the  wind  axes  and  the  vehicle  body  axes  of  a 
rigid  body  defines  the  angle  of  attack,  a,  and  the  sideslip  angle,  0.  These 
angles  are  convenient  independent  variables  for  use  in  the  expression  of 
aerodynamic  force  and  moment  coefficients. 

Wind  axes  are  not  generally  used  in  the  analysis  of  the  motion  of  a  rigid 
body,  because,  as  in  the  case  of  the  earth  axes,  the  moment  of  inertia  and 
product  of  inertia  terms  in  the  three  rotational  equations  of  motion  vary  with 
time,  a,  and  $. 

4.5  DERIVATION  OF  THE  RIGHT  HAND  SIDE  (RHS)  OF  THE  EQUATIONS  OF  MOTION 

The  RHS  of  the  equation  represents  the  aircraft  response  to  forces  or 
moments.  Through  the  application  of  Newton's  second  law,  two  vector  relations 
can  be  used  to  derive  the  six  required  equations,  three  translational  and 
three  rotational.  The  actual  aircraft  will  be  flexible,  which  gives  rise  to 
aeroelastic  effects  and  are  additional  degrees  of  freedom  (requiring  more 
variables  and  equations).  These  will  be  considered  separately  (in  Chapter  12). 
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In  order  to  derive  the  equations  of  motion  we  will  make  the  following 

ASSUMPTION;  The  aircraft  is  a  rigid  body. 

As  noted  earlier,  Newton's  second  law  is  valid  only  in  an  inertial 
coordinate  system.  For  most  aircraft,  the  fixed  earth  axis  system  can  be 
assumed  to  be  an  inertial  coordinate  system.  In  order  to  do  this,  we  make  the 
following 

ASSUMPTION:  The  earth  and  atmosphere  are  fixed  in  inertial  space. 

In  addition,  most  motion  of  interest  in  stability  and  control  takes  place 
in  a  relatively  short  time.  We  can  also,  therefore,  usually  make  the  following 

ASSUMPTION:  Mass  (m)  is  constant  (dm/dt  =  0). 

4.5.1  Translational  Force  Relations 

The  vector  equation  for  the  aircraft  translation  from  Newton's  second  law 
(eq.  4.1)  is 
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=  d(mVT) 


dt 


XYZ 


(4.9) 


where  VT  is  the  true  velocity  of  the  aircraft.  Figure  4.6  shows  how  this 
vector  changes  in  both  magnitude  and  direction  with  respect  to  the  xyz  (body) 
and  XYZ  (fixed  earth)  axes. 


From  vector  analysis,  the  derivative  of  the  velocity  VT  in  the  inertial 
(fixed  earth)  coordinate  system  is  related  to  the  derivative  of  VT  in  the  body 
axis  system  through  the  relationship 
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FIGURE  4.6.  TRUE  VELOCITY  IN  BODY  AND  FIXED  EARTH  AXES 
COORDINATE  SYSTEMS 


Substituting  this  into  Equation  4.9  (and  assuming  mass  is  constant),  the 
applied  force  is 
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VT  and  to  are  two  of  the  four  vectors  used  in  the  equations  of  motion  to 
describe  the  vehicle  motion  (F  and  G  are  the  other  two).  They  are  defined  as 
follows: 


VT  =  Ul  +  V3  +  w£  (4.12) 

where 

U  =  forward  velocity 
V  =  side  velocity 
W  =  vertical  velocity 


and 


a)  *=  Pi  +  Qj  +  RJc 


(4.13) 


where 

P  •  roll  rate 
Q  -  pitch  rate 
R  «  yaw  rate 

The  relationship  of  the  true  velocity  and  its  components  to  o  and  3  and 
the  body  axis  coordinate  system  is  shown  in  Figure  H.l. 
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FIGURE  4.7.  VELOCITY  COMPONENTS  AND  THE  AERODYNAMIC 
ORIENTATION  ANGLES,  a  AND  3 


The  angles  a  and  3  can  be  expressed  in  terms  of  the  velocity  components  as 
follows: 


Sin  a  =  - - 

VTcos  3 

If  3  is  small  (ASSUMPTION) ,  then  cos  3  «  1  and 


(4.14) 


sin  a  = 

VT 

If  a  is  also  small  (ASSUMPTION) ,  then 
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For  angle  of  sideslip 
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If  0  is  small,  then  sin  0=0  and 

0  « 


v 


V 

V„ 


(4.17) 


(4.18) 


Using  equations  4.12  and  4.13,  the  translational  equation  (4.11)  can  now 
be  written  in  component  form  as 

;  i  5  k 

F  =  m  [  Ui  +  vj  +  Wk  + 
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Expanding 

F  =  m  [Ui  +  Vj  +  Wk  +  (QW  -  RV)T  -  (PW  -  RU)]  +  (PV  -  QU)Jt] 
Rearranging 

F  =  m  [(U  +  QW  -  RV)I  +  (V  +  RU  -  PW)]  +  (W  +  PV  -  QU)k] 

In  component  form,  the  sum  of  forces  in  the  body  axis  system  is 

F  =  fI  +  fi  +  fE 

x  y  J  z 

this  results  in  three  component  translational  equations: 
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Fy  =  m  (V  +  RU  -  PW) 


Fz  =  m  (W  +  PV  -  QU) 


(4.23) 

(4.24) 

(4.25) 


4.5.2  Rotational  Equations 

Once  again  from  Newton's  second  law  (eq.  4.2), 
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Equation  4.26  states  the  change  in  angular  momentum,  H,  is  equal  to  the 
applied  moments,  G. 

Angular  momentum  should  not  be  as  difficult  to  understand  as  some  people 
would  like  to  make  it.  It  can  be  thought  of  as  linear  momentum  with  a  moment 
arm  included.  Consider  a  ball  swinging  on  the  end  of  a  string,  at  any  instant 
of  time,  as  shown  in  Figure  4.8. 
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FIGURE  4.8.  ANGULAR  MOMENTUM 
The  linear  momentum  of  this  system  would  be: 

Linear  Momentum  =  mV 

Angular  momentum  is  defined  as  H,  where  H  =  r  X  Linear  Momentum  and, 

since  in  the  example  of  Figure  4.8,  the  angle  between  r  and  V  is  90 
degrees,  the  magnitude  of  the  angular  momentum  is  mrV. 

Just  as  a  force  F  changes  linear  momentum,  (F  =  ^  mV) ,  a  moment  G  will 

__  -3 

change  angular  momentum  (G  =  ^  H).  A  moment  is  related  to  a  force  in  the 
same  manner  that  angular  momentum  is  related  to  linear  momentum: 

Moment  =  r  X  Force 

Angular  Momentum  =  r  X  Linear  Momentum 

In  order  for  us  to  determine  the  angular  momentum  of  the  aircraft, 
consider  a  small  element  of  mass  ,  somewhere  in  the  aircraft,  a 

distance  r1  from  the  eg  (Figure  4.9). 


4.13 


FIGURE  4.9.  ELEMENTAL  DEVELOPMENT  OF  RIGID 
BODY  ANGULAR  MOMENTUM 

The  angular  momentum  of  is 

K  -  ^  XmV,  «  (?,  X  VJ  (4.27) 

and 

( i.e. ,  in  the  inertial  coordinate  system) 

XYZ 

Again  from  vector  analysis,  the  rate  of  change  of  the  radius  vector  r  can  be 
related  to  the  body  axis  system  (xyz)  by 


since  the  aircraft  is  a  rigid  body  rx  does  not  change  with  time  (assuming  no 
aeroelastic  effects).  Therefore,  the  first  term  can  be  excluded,  and  the 
inertial  velocity  of  the  element  rr^  is 
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i  5  k 
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( Qz-Ry )  ( Rx-Pz )  ( Py-Qx ) 

So  the  components  of  H  are 

H  =  Jv  PAy(Py  -  Qx)dv  -  Xv  paz(Rx  -  Pz)dV  (4.37) 

Hy  =  K  paz(Qz  -  Ry)dV  -  Jv  pAx(Py  -  Qx)dV  (4.38) 

Hz  =  K  pax(Rx  “  "  Xv  PAy(Qz  -  Ry)dV  (4.39) 

Rearranging  the  equations 

H  =  p  lv  PA(y2  +  z2  )dV  -  Q  Xv  PA xydv  -  R  Jy  pAxzdV  (4.40) 
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Hy  -  Q  Jv  pA(x2  +  z2)dV  -  R  Jv  pAyzdV  -  P  Xv  PAxydV  (4.41) 

H*  “  R  Xv  pa(x2  +  y2  )ctV  -  P  Jv  pAxzdV  -  Q  Jv  pAyzdV  (4.42) 

The  integrals  are  now  recognizable  as  moments  and  products  of  inertia.  The 
moments  of  inertia  are  defined  as 
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-  Jv 

PA(X2  +  z2  )dv 

(4.44) 

■  Jv 

pA(x2  +  y2  )dV 

(4.45) 

These  are  a  measure  of  resistance 

to  rotation  and  are  never  zero. 

They  are 

illustrated  in  Figure  4.10. 


FIGURE  4.10.  MOMENT  OF  INERTIA  (Iy) 
The  products  of  inertia  are  defined  as  (Figure  4.11) 


*xy 

=  *yx 

=  Jv 

PAxy  dv 

(4.46) 

=  Jxy 

-  Jv 

pAyz  dv 

(4.47) 

- 

-  Jv 

pAxz  dv 

(4.48) 
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FIGURE  4.11.  PRODUCT  OF  INERTIA  (Ixx ) 

Products  of  inertia  are  measures  of  symmetry.  They  are  zero  for  views  having 
a  plane  of  symmetry. 

Substituting  into  equations  4.40  to  4.42,  we  find  the  angular  momentum  of 
a  rigid  body  is 


So  that 


H  ■=  hI  +  h-T  +  hJc 

x  y J  z 


H 

«=  PI 

-  QI 

-  RI 

X 

X 

~  xy 

xz 

H 

«  QI 

-  RI 

-  PI 

y 

y 

yz 

xy 

H 

«  RI 

-  PI 

-  QI 

z 

z 

x  z 

yz 

(4.49) 

(4.50) 

(4.51) 

(4.52) 


An  aircraft  is  normally  symmetric  about  the  xz-axis  as  illustrated  in 
Figure  4.12.  In  order  to  simplify  the  RHS  of  the  equations  of  motion, 
therefore,  we  normally  make  the  following 


ASSUMPTION;  Hie  xz-plane  is  a  plane  of  symmetry. 

This  causes  two  products  of  inertia,  I  and  I  to  be  zero.  These  may  be 
cancelled  out  of  the  equations  of  motion.  This  restriction  of  the  equations 
of  motion  (xz-plane  symmetry)  can  be  easily  removed  by  including  these  terms. 
With  the  assumption,  the  angular  momentum  of  a  symmetric  aircraft  simplifies  to 

H  =  (PIx  -  RIxz)  I  +  QIy3  +  (RIZ  -  PIX2)  E  (4.53) 
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*  * 

FIGURE  4.12.  AIRCRAFT  INERTIAL  PROPERTIES  WITH 
AN  x-z  PLANE  OF  SYMMETRY 

The  equation  for  angular  momentum  can  now  be  substituted  into  the  moment 
equation.  Remember 
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Since  the  body  axis  system  is  used,  the  moments  of  inertia  and  the 
products  of  inertia  are  constant.  Therefore,  by  differentiating  and 
substituting,  the  moment  equation  becomes 


G  =  (PI,  -  RI„>  i  +  QIy j  *  (FI,  -  fl„)  E  + 


1 

p 


j 

Q 


E 

R 


<PIx-RIxZ>  QIy  (RIZ-PIXZ) 


Therefore,  the  rotational 

component  equations  are. 

(4.58) 

G 

X 

=  PI 

X 

+  QR  (I  -  I  )  - 

z  y 

(R  +  PQ)  I 

xz 

(4.59) 

Gy 

=  Ql 

y 

-  PR  (I  -  I  )  + 

(P2  -R2)  Ixz 

(4.60) 

G 

Z 

=  RI 

z 

+  PQ  (Iy  -  Ix)  + 

(QR  -  P)  Ixz 

(4.61) 

This  completes  the  development  of  the  RHS  of  the  six  equations  (equations  4.23 
to  4.25,  and  4.59  to  4.61). 

4.6  DERIVATION  OF  THE  LHS  OF  THE  EQUATIONS  OF  MOTION 

The  equations  of  motion  relate  the  vehicle  motion  to  the  applied  forces 
and  moments: 


LHS  RHS 

Applied  Forces  and  Moments  =  Observed  Vehicle  Motion 

Fx  =  m(U  +  QW  -  FV) 

Gx  =  K  +  "  V  "  <R  +  PQ)lyz 

etc. 

The  RHS  of  each  of  these  six  equations  has  been  completely  expanded  in  terms 
of  easily  measured  quantities.  The  LHS  must  also  be  expanded  in  terms  of 
convenient  variables.  In  order  to  do  this,  we  must  be  able  to  relate  the 
orientation  of  the  body  axes  (xyz)  to  the  moving  earth  axes  (XYZ).  This  is 


done  through  the  use  of  Euler  angles.  The  moving  earth  axis  system  is  used 
because  we  will  be  concerned  with  the  orientation  of  the  aircraft  with  respect 
to  the  earth  and  not  its  position  (location  of  the  eg)  with  respect  to  the 
earth. 


4.6.1  Euler  Angles 

The  orientation  of  any  coordinate  system  relative  to  another  can  be  given 
by  three  angles  (Euler  angles),  which  are  consecutive  rotations  about  the  z, 
y,  and  x  axes,  in  that  order,  that  carry  one  frame  into  coincidence  with  the 
other.  In  flight  dynamics,  the  Euler  angles  used  are  those  which  rotate  the 
vehicle  carried  moving  earth  axis  system  into  coincidence  with  the  relevant 


vehicle  axis  system  (Figure  4.13). 


FIGURE  4.13.  THE  EULER  ANGLE  ROTATIONS 
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Euler  angles  are  expressed  as  YAW  (vd,  PITCH  (0),  and  ROLL  The 
sequence  (YAW,  PITCH,  ROLL)  must  be  maintained  to  arrive  at  the  proper  set  of 
Euler  angles.  The  Euler  angles  are  defined  as  follows: 


V  -  Yaw  Angle  - 


The  angle  between  the  projection  of  x  vehicle  axes 
onto  the  horizontal  plane  and  the  initial  reference 
position  of  the  X  earth  axis.  (Yaw  angle  is  the 
vehicle  heading  only  if  the  initial  reference  is 
North) . 


0  -  Pitch  Angle  -  The  angle  measured  in  a  vertical  plane  between  the  x 

vehicle  axis  and  the  horizontal  plane. 


<(>  ~  Roll  Angle  -  The  angle,  measured  in  the  yz  plane  of  the  vehicle 

system,  between  the  y  axis  and  the  horizontal  plane. 
This  is  the  same  as  bank  angle  for  a  given  and  0, 
and  is  a  measure  of  the  rotation  about  the  x  axis  to 
put  the  aircraft  in  the  desired  position  from  a 
wing's  horizontal  condition. 

The  accepted  limits  on  the  Euler  angles  are: 

-180°  <  <|>  <  +  180° 

-90°  <  0  <  +  90° 

-180°  <  <f>  <  +  180° 

The  importance  of  the  sequence  of  the  Euler  angle  rotations  cannot  be 
overemphasized.  Finite  angular  displacements  do  not  behave  as  vectors. 
Therefore,  if  the  sequence  is  performed  in  a  different  order  than  \f/,  0,  <p,  the 
final  result  will  be  different.  This  fact  is  clearly  illustrated  by  the  final 
aircraft  attitudes  in  Figure  4.14  in  which  two  rotations  of  equal  magnitude 
have  been  performed  about  the  x  and  y  axes,  but  in  opposite  order.  Addition 
of  a  rotation  about  a  third  axis  does  nothing  to  improve  the  outcome . 

Euler  angles  are  very  useful  in  describing  the  orientation  of  flight 
vehicles  with  respect  to  inertial  space.  Consequently,  angular  rates  in  an 
inertial  system  ($,  ©,  <J/)  can  be  transformed  to  angular  rates  in  the  vehicle 
axes  (P,  Q,  R)  using  Euler  angle  transformations  as  developed  in  the  next 
subsection. 
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ROTATION  SEQUENCE  1 


ROTATION  SEQUENCE  2 


FIGURE  4.14.  DEMONSTRATION  THAT  FINITE  ANGULAR 

DISPLACEMENTS  DO  NOT  BEHAVE  AS  VECTORS 


4.6.2  Angular  Velocity  Transformations 

We  need  to  develop  equations  to  transform  the  angular  rates  from  the 

moving  earth  axis  system  (^,  0,  ♦)  into  angular  rates  about  the  vehicle  axis 

system  (P,  Q,  R)  for  any  aircraft  attitude.  The  derivation  (by  vector 

resolution)  is  presented  in  the  following  paragraphs. 

It  is  easy  to  see  that  when  an  aircraft  is  pitched  up  and  banked,  the 

vector  ip  will  have  components  along  the  x,  y,  and  z  body  axes  (Figure  4.15). 

Remember,  is  the  angular  velocity  about  the  Z  axis  of  the  moving  earth  axis 

system  (it  can  be  thought  of  as  the  rate  of  change  of  aircraft  heading). 

Although  it  is  not  shown  in  Figure  4.15,  the  aircraft  may  have  a  value  of  0 

and  $.  In  order  to  derive  the  transformation  equations,  it  is  easier  to 

analyze  one  vector  at  a  time.  First  resolve  the  components  of  ip  on  the  body 

•  • 

axes.  Then  do  the  same  with  0  and  i.  The  components  can  then  be  added  and 
the  total  transformation  will  result. 


FIGURE  4.15.  COMPONENTS  OF  ip  ALONG  x,  y,  AND  z  BODY  AXES. 

(NOTE:  THE  X  AND  Y  AXES  OF  THE  MOVING  EARTH 
AXIS  SYSTEM  ARE  NOT  SHOWN. ) 


Step  1  -  Resolve  the  components  of  ip  along  the  body  axes  for  any 
aircraft  attitude. 

It  is  easy  to  see  how  ip  reflects  to  the  body  axis  by  starting  with 
an  aircraft  in  straight  and  level  flight  and  changing  the  aircraft  attitude 
one  angle  at  a  time.  In  keeping  with  convention,  the  sequence  of  change  will 
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be  yaw,  pitch  and  bank. 

First,  it  can  be  seen  from  Figure  4.16  that  the  Z-axis  of  the  moving 
earth  axis  system  remains  aligned  with  the  z-axis  of  the  body  axis  system 
regardless  of  the  angle  ^  if  0  and  $  are  zero;  therefore,  \j/  does  not  affect  P, 
and  Q. 

.*.  R  =  (when  0  =  <(>  =  0) 


FIGURE  4.16.  DEVELOPMENT  OF  AIRCRAFT  ANGULAR  VELOCITIES 
BY  THE  EULER  ANGLE  YAW  RATE  (\p  ROTATION) 

Next,  consider  pitch  up.  In  this  attitude,  ij)  has  components  on  the  x  and 
z-body  axes  as  shown  in  Figure  4.17.  As  a  result,  ^  will  contribute  to  the 
angular  rates  about  these  axis. 

P  =  -  \p  sin  ©  (4.62) 

R  =  cos  0  (4.63) 

With  just  pitch,  the  Z-axis  remains  perpendicular  to  the  y-body  axis,  so  Q  is 
not  affected  by  in  this  attitude. 

Next,  bank  the  aircraft,  leaving  the  pitch  as  it  is  (Figure  4.18). 


FIGURE  4.17.  DEVELOPMENT  OF  AIRCRAFT  ANGULAR  VELOCITIES 
BY  THE  EULER  ANGLE  YAW  RATE  (6  ROTATION) 


FIGURE  4.18.  DEVELOPMENT  OF  AIRCRAFT  ANGULAR  VELOCITIES 
BY  THE  EULER  ANGLE  YAW  RATE  (♦  ROTATION) 
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All  of  the  components  are  now  illustrated.  Notice  that  roll  did  not 
change  the  effect  of  *  on  P.  Hie  components,  therefore,  of  in  the  body 
axes  for  any  aircraft  attitude  are 

p  «  -  ^  sin  0  (4.64) 

(Effect  of  ij»  only)  Q  *  i  c°s  ©  s:i-n  ♦  (4.65) 

r  m  $  cos  0  cos  4>  (4.66) 

Step  2  -  Resolve  the  components  of  6  along  the  body  axes  for  any 
aircraft  attitude. 

Remember,  0  is  the  angle  between  the  x-body  axis  and  the  local  horizontal 
(Figure  4.19).  Once  again,  change  the  aircraft  attitude^  by  steps  in  the 
sequence  of  yaw,  pitch,  and  bank  and  analyze  the  effects  of  0. 


z 

FIGURE  4.19.  CONTRIBUTION  OF  THE  EULER  PITCH  ANGLE  RATE 
TO  AIRCRAFT  ANGULAR  VELOCITIES  (©  ROTATION) 

It  can  be  seen  immediately  that  the  yaw  angle  has  no  effect.  Likewise 
when  pitched  up,  the  y-body  axis  remains  in  the  horizontal  plane.  Therefore, 
§  is  the  same  as  Q  in  this  attitude  and  the  component  is  equal  to 

Q  =  0 

Now  bank  the  aircraft. 
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FIGURE  4.20.  CONTRIBUTION  OF  THE  EULER  PITCH  ANGLE  RATE 
TO  AIRCRAFT  ANGULAR  VELOCITIES  ( <j>  ROTATION) 

It  can  be  seen  from  Figure  4.20  that  the  components  of  9  on  the  body  axes  are 


Q  =  9  cos  <j>  (4.67) 

R  =  -9  sin  <p  (4.68) 

Notice  that  P  is  not  affected  by  9  since  by  definition  9  is  measured  on  an 
axis  perpendicular  to  the  x  body  axis. 

steP  3  ~  Resolve  the  components  of  $  along  the  body  axes. 

This  one  is  easy  since  by  definition  is  measured  along  the  x  body  axis. 
Therefore,  i  affects  the  value  of  P  only,  or 

p  =  ♦  (4.69) 


•  •  • 

The  components  of  ip,  9,  and  <f>  along  the  x,  y,  and  z  body  axes  for  any 
^i^craft  attitude  have  been  derived.  These  can  now  be  summed  to  give  the 
transformation  equations. 
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P  «  sin  0)1 


(4.70) 


Q  «  (0  cos  ♦  +  'i'  cos  0  sin  4>)] 


(4.71) 


R  «  (if/  cos  0  cos  ♦  -  0  sin  $>)£  (4.72) 

With  these  equations  it  is  now  possible  to  transform  the  equations  of 
motion  written  in  bocty  axis  terms  (U,  V,  W,  P,  Q,  and  R)  in  terms  of  the 
motion  seen  in  the  inertial  (earth  axis)  system  (u,  V,  W,  0/  and  4>).  In 
that  case,  the  resulting  equations  are  six  simultaneous  nonlinear  differential 
equations  that  are  of  the  first  order  in  U,  V,  and  W  and  of  the  second  order 
in  ij>,  0,  and  $.  In  order  to  completely  describe  the  trajectory  of  the 
aircraft  in  the  inertial  coordinate  system,  a  similar  transformation  is 
required  to  relate  U,  V,  and  W  to  the  velocities  in  the  inertial  coordinate 
system  (not  covered  in  this  text). 

Equations  4.70  to  4.72  are  known  as  the  parametric  equations  and  they, 
along  with  the  six  equations  of  motion,  can  be  used  to  describe  the  complete 
motion  of  the  aircraft. 


4.6.3  Initial  Breakdown  of  the  LHS 

In  general,  the  applied  forces  and  moments  on  the  LHS  can  be  broken  up 
according  to  the  sources  shown  below. 


SOURCE 

Aero¬ 

dynamic 

Direct 

Thrust 

Gravity 

Gyro¬ 

scopic 

Other 
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xoth 
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^oth 
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mt 

0 

flyro 

it 

HI  h 

O 

Ya 

Y„ 

0 

Yoth 

Lt 

0 

^"oth 

□  CC 

5 

Na 

nt 

0 

N„ 

gyro 

=  mil H -  (4.73) 

=  mW  +  —  (4.73a) 
=  Qly  +  -  -  -  (4.74) 

=  mV  +  *--  (4.74a) 
=  P'.  +  -"  (4.75) 
=  *‘«+  (4.75a) 
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1.  Aerodynamic  Forces  and  Moments  -  These  will  be  further  expanded  into 
stability  parameters  and  derivatives  (discussed  in  subsection  4.6.4). 

2.  Direct  Thrust  Forces  and  Moments  -  These  terms  include  the  effect  of 
the  thrust  vector  itself  -  they  usually  do  not  include  the  indirect 
or  induced  effects  of  jet  flow  or  running  propellers  (discussed  in 
subsection  4.6.5). 

3.  Gravity  Forces  —  These  vary  with  orientation  of  the  gravity  vector 
(discussed  in  subsection  4.6.6). 

4.  Gyroscopic  Moments  -  These  occur  as  a  result  of  large  rotating 
masses  such  as  engines  and  props  (discussed  in  section  4.6.7). 

5.  Other  Sources  -  These  include  spin  chutes,  reaction  controls,  etc. 
(not  discussed  in  this  chapter). 

4.6.4  Aerodynamic  Forces  And  Moments 

By  far  the  most  important  forces  and  moments  on  the  LHS  of  the  equation 
are  the  aerodynamic  terms.  Unfortunately,  they  are  also  the  most  complex.  As 
a  result,  certain  simplifying  assumptions  are  made,  and  several  of  the  smaller 
terms  are  arbitrarily  excluded  to  simplify  the  analysis.  Remember  we  are  not 
trying  to  design  an  aircraft  around  some  critical  criteria.  We  are  only 
trying  to  derive  a  set  of  equations  that  will  help  us  analyze  the  important 
factors  affecting  aircraft  stability  and  control. 

4. 6. 4.1  Choice  Of  Axis  System.  Consider  only  the  aerodynamic  forces  on  an 
aircraft.  Summing  forces  along  the  x  body  axis  (Figure  4.21) 

Fx  «*  L  sin  a  -  D  cos  a  (4.76) 


FIGURE  4.21.  CHOICE  OF  AXIS  SYSTEM 
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Notice  that  if  the  forces  were  summed  along  the  xs  stability  axis  (Figure 
4.21) ,  it  would  be 


Fx  =  -  D  (4.77) 

It  would  simplify  things  if  the  stability  axes  were  used  for  development  of 
the  aerodynamic  forces.  A  small  angle  assumption  enables  us  to  do  this: 

cos  a  «  1 

sin  a  «  0 

Using  this  assumption,  equation  4.76  reduces  to  equation  4.77.  Whether  it  be 
thought  of  as  a  small  angle  assumption  or  as  an  arbitrary  choice  of  the 
stability  axis  system,  the  result  is  less  complexity.  This  would  not  be  done 
for  preliminary  design  analyses;  however,  for  the  purpose  of  deriving  a  set  of 
equations  to  be  used  as  an  analytical  tool  in  determining  handling  qualities, 
the  assumption  is  perfectly  valid,  and  is  surprisingly  accurate  for  relatively 
large  values  of  a.  It  should  be  noted  that  lift  and  drag  are  defined  to  be 
positive  as  illustrated.  Thus  these  quantities  have  a  negative  sense  with 
respect  to  the  stability  axis  system. 

The  aerodynamic  terms  will  be  developed  using  the  stability  axis  system 


so  that  the  equations  assume  the  form, 

"DRAG"  -D  +  5^  +  Xg  +  Xoth  =  mU  + -  (4.78) 

"LIFT"  -L  +  ZT  +  Zg  +  Zoth  =  mW  + -  (4.79) 

"PITCH"  M.  +  ML  +  M  +  M  ..  =  QI  + -  (4.80) 

"SIDE"  Ya  +  Yt  +  Yg  +  Yoth  =  mV  + -  (4.81) 

"ROLL"  L+L+L  +  L  v>  =  PI  +  -  (4.82) 

a  gyro  oth  x  '  ' 

"YAW"  N  +  N  +  N  +  N  .  .  =  RI  + -  (4.83) 

A  T  gy  r  o  o  t  h  z 


4.6.4.1A  Coordinate  Systems  and  Transformations,  The  five  orthogonal  coordinate 
systems  are  related  by  the  figure  below.  In  general,  you  need  only  to  rotate 
through  the  given  angles  to  transfer  from  one  set  of  axes  to  another. 
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FIGURE  4.21A  AXIS  SYSTEM  RELATIONSHIPS 


It  is  often  convenient  to  measure  forces  or  moments  in  a  certain  axis  system  while 
the  equations  of  motion  are  better  understood  in  another.  If  you  are  transforming 
in  the  direction  indicated,  use  the  normal  rotation  matrices.  If  transforming 
opposite  to  the  indicated  direction  use  the  transpose  of  the  rotation  matrices.  In 
either  case  put  in  the  measured  angle(s)  without  changing  sign(s) .  For  example, 
weight  is  measured  easily  in  the  earth  axis  system.  To  transform  to  the  body  axis 
the  process  is  to  pre-multiply  the  weight  vector (earth  axis)  by  R^  then  R^  and 
then  R^  as  shown. 


1 


Weight 


J  Body 


=  R: (0A)  R2(S2)  R 3(Y3) 


0 

0 


Weight 


Earth 


Where  the  rotation  matrices  are 


1  0 

0 
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0 
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0 
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0 

1 

0 
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0 

X 
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0 
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■J 

0 

0 

1 

And  the  (*)  represents  the  appropriate  angle  from  figure  4.21A 


4. 6. 4. 2  Expansion  of  Aerodynamic  Terms.  A  stability  and  control  analysis  is 
concerned  with  how  a  vehicle  responds  to  perturbation  inputs.  For  instance, 
up  elevator  should  cause  the  nose  to  come  up;  or  for  turbulence  caused 
sideslip,  the  aircraft  should  realign  itself  with  the  relative  wind. 
Intuitively,  the  aerodynamic  terms  have  the  most  effect  on  the  resulting 
motion  of  the  aircraft.  Unfortunately,  the  above  equations  (that  result  from 
summing  forces  and  moments),  are  non-linear,  and  exact  solutions  are 
impossible.  In  view  of  the  complexity  of  the  problem,  linearization  of  the 
equations  brings  about  especially  desirable  simplifications.  The  linearized 
model  is  based  on  the  assumption  of  small  disturbances  and  the  small 
perturbation  theory.  This  model,  nonetheless,  gives  quite  adequate  results 
for  engineering  purposes  over  a  wide  range  of  applications;  because  the  major 
aerodynamic  effects  are  nearly  linear  functions  of  the  variables  of  interest, 
and  because  quite  large  disturbances  in  flight  may  correspond  to  relatively 
small  disturbances  in  the  linear  and  angular  velocities. 

4. 6. 4. 3  Small  Perturbation  Theory.  The  small  perturbation  theory  is  based  on 
a  simple  technique  used  for  linearizing  a  set  of  differential  equations.  In 
aircraft  flight  dynamics,  the  aerodynamic  forces  and  moments  are  assumed  to  be 
functions  of  the  instantaneous  values  of  the  perturbation  velocities,  control 
deflections,  and  of  their  derivatives.  They  are  obtained  in  the  form  of  a 
Taylor  series  in  these  variables,  and  the  expressions  are  linearized  by 
excluding  all  higher-order  terms.  To  fully  understand  the  derivation,  some 
assumptions  and  definitions  must  first  be  established. 

4. 6. 4. 3.1  The  Small  Disturbance  Assumption  -  A  summary  of  the  major 
variables  that  affect  the  aerodynamic  characteristics  of  a  rigid  body  or  a 
vehicle  is  given  below. 

1.  Velocity,  temperature,  and  altitude:  These  variables  may  be 
considered  directly  or  indirectly  functions  of  Mach,  Reynolds 
number,  and  dynamic  pressure.  Velocity  may  be  resolved  into 
components  U,  V,  and  W  along  the  vehicle  body  axes. 

2.  Angle  of  attack,  a,  and  angle  of  sideslip,  (3:  These  variables 
may  be  used  with  the  magnitude  of  the  total  velocity,  V  ,  to 
express  the  orthogonal  velocity  components  U,  V,  and  W.  TIt  is 
more  convenient  to  express  variation  of  force  and  moment 
characteristics  with  these  angles  as  independent  variables 
rather  than  the  velocity  components. 
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3.  Angular  velocity:  This  is  usually  resolved  into  components,  P, 
Q,  and  R  about  the  vehicle  body  axes. 

4.  Control  surface  deflections:  These  are  used  primarily  to 

change  or  balance  aerodynamic  forces  and  moments,  and  are 
accounted  for  by  S  ,  $  ,  5r  (the  elevator,  aileron  and  rudder 
deflection,  respectively) . 

Because  air  has  mass,  the  flow  field  cannot  adjust  instantaneously  to  sudden 
changes  in  these  variables,  and  transient  conditions  exist.  In  some  cases, 
these  transient  effects  become  significant.  Analysis  of  certain  unsteady 
motions  may  therefore  require  consideration  of  the  time  derivatives  of  the 
variables  listed  above.  In  other  words: 


VARIABLE 

FIRST  DERIVATIVE 

SECOND  DERIVATIVE 

D 

U  a  P 

U  i  (3 

**  ,,  »• 

U  a  (3 

L 

Ma 

Are  a 

«  Function  < 

P  Q  R 

•  •  • 

P  Q  R 

P  Q  R 

Ya 

L 

of 

Se  S  5 

ear 

i  1  5 

ear 

•  '  *  *  it 

S  5  S 

ear 

A 

Na  - 

p  M  Re  T 

assumed  constant 

This  rather  formidable  list  can  be  reduced  to  workable  proportions  by 
assuming  that  the  vehicle  motion  consists  only  of  small  deviations  from  some 
initial  reference  condition.  In  addition,  a  Taylor  series  expansion,  with 
higher  order  terms  assumed  negligible,  is  used  to  determine  the  effect  of 
these  small  perturbations  on  the  aircraft.  Fortunately,  this  small 
disturbance  assumption  applies  to  many  cases  of  practical  interest  and,  as  a 
bonus,  stability  parameters  and  derivatives  derived  under  this  assumption 
continue  to  give  good  results  for  somewhat  larger  motions. 

The  variables  are  considered  to  consist  of  some  equilibrium  value  plus  an 
incremental  change,  called  the  "perturbed  value."  The  notation  for  these 
perturbed  values  is  usually  lower  case.  For  example, 

P  =  P0  +  p 

U  =  U0  +  u 
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In  summary,  the  small  disturbance  assumption  is  applied  in  three  steps: 
assuming  an  initial  (equilibrium)  condition  (described  in  subsection 
4. 6. 4. 3. 3),  assuming  vehicle  motion  consists  of  small  perturbations  about  this 
condition  and,  using  a  first  order  Taylor  series  expansion  (described  in 
subsection  4. 6. 4. 3. 4)  to  determine  the  effect  of  these  small  perturbations. 
As  an  additional  consequence,  the  small  perturbation  assumption  allows  us  to 
decouple  the  longitudinal  and  lateral-directional  equations  as  discussed  in 
the  next  subsection. 

4. 6. 4. 3. 2  Longitudinal  and  Lateral-Directional  Equations  -  It  has  been 
found  from  experience  that,  when  operating  under  the  small  perturbation 
assumption,  the  vehicle  motion  can  be  thought  of  as  two  independent 
(decoupled)  motions,  each  of  which  is  a  function  only  of  the  variables  shown 
below. 


1.  Longitudinal  Motion 

(D,  L,  Ma)  =  f  (U,  a,  k,  Q,  Se)  (4.84) 

2.  Lateral-Directional  Motion 

(Ya,  La,  Na)  =  f  (8,  8,  p,  R,  $a,  sr)  (4.85) 

The  equations  are  grouped  and  named  in  the  above  manner  because  the  state 
variables  of  the  first  group  U,  a,  a,  Q,  Se  are  known  as  the  longitudinal 

variables  and  those  of  the  second  group,  8/  0,  P,  R,  Sa ,  and  &r ,  are  known 

as  the  lateral-directional  variables.  With  the  conventional  simplifying 
assumptions,  the  longitudinal  and  lateral-directional  variables  will  appear 
explicitly  only  in  their  respective  group.  This  separation  will  also  be 

displayed  in  the  aerodynamic  force  and  moment  terms  and  the  equations  will 
completely  decouple  into  two  independent  sets. 

4. 6. 4. 3. 3  Initial  Conditions  -  As  stated  earlier,  we  will  assume  that 
the  motion  consists  of  small  perturbations  about  some  initial  equilibrium 
condition.  The  condition  we  will  assume  is  steady  straight  symmetrical 
flight.  This  condition  is  a  combination  of  the  following  motions: 

Steady  Flight.  Motion  with  zero  rates  of  change  of  the  linear  and 
angular  velocity  components,  i.e., 

U  =  V  =  W=  P  =  Q  =  R  =  0. 
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Straight  Flight.  Motion  with  zero  angular  velocity  components ,  P,  Q,  and 

R  =  0. 

Symmetric  Flight.  Motion  in  which  the  vehicle  plane  of  symmetry  remains 
fixed  in  space  throughout  the  maneuver.  The  unsymmetric  variables  P,  R,  V,  <f>, 
and  (3  are  all  zero  in  symmetric  flight.  Some  symmetric  flight  conditions  are 
wings-level  dives,  climbs,  and  pull-ups  with  no  sideslip.  In  steady  straight 
symmetric  flight,  the  aircraft  is  assumed  to  be  flying  wings  level  with  all 
components  of  velocity  zero  except  U0  and  w0 .  Therefore,  with  reference  to 
the  body  axis 

VT  =  UQ  *  constant 

WQ  =  small  constant  .*.  =  small  constant 

V0  *  0  :.  e0  =0 

po  “  Q0  “  R0  “  0 

We  have  already  found  that  the  equations  of  motion  simplify  considerably 
when  the  stability  axis  is  used  as  the  reference  axis.  This  idea  will  again 
be  employed  and  the  final  set  of  boundary  conditions  will  result.  This, 
therefore,  is  another 

ASSUMPTION:  VT  =  UQ  =  constant 

Wo  “  0  «d  *  0 

V0  -  o  .-.  e0  -  0 

Po  -  Qo  -  Ro  “  0 

(p,  M,  Re,  aircraft  configuration)  =  constant 

4. 6. 4. 3. 4  Expansion  By  Taylor  Series.  As  stated  earlier,  the  equations 
resulting  from  summing  forces  and  moments  are  nonlinear  and  exact  solutions 
are  not  obtainable.  An  approximate  solution  is  found  by  linearizing  these 
equations  using  a  Taylor  Series  expansion  and  neglecting  higher  ordered  terms. 
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As  an  introduction  to  this  technique,  assume  some  arbitrary  non-linear 
function,  f(U),  having  the  graphical  representation  shown  in  Figure  4.22. 


FIGURE  4.22.  APPROXIMATION  OF  AN  ARBITRARY 
FUNCTION  BY  TAYLOR  SERIES 

A  Taylor  Series  expansion  will  approximate  the  curve  over  a  short  span,  AU. 
Ihe  first  derivative  assumes  the  function  between  AU  to  be  a  straight  line 
with  slope  3f (U0  )/3U.  This  approximation  is  illustrated  in  Figure  4.23. 


FIGURE  4.23.  FIRST  ORDER  APPROXIMATION  BY  TAYLOR  SERIES 
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To  refine  the  accuracy  of  the  approximation,  a  second  derivative  term  is 
added,  The  second  order  approximation  is  shown  in  Figure  4.24. 


Uc+  AU 


FIGURE  4.24.  SECOND  ORDER  APPROXIMATION  BY  TAYLOR  SERIES 

Additional  accuracy  can  be  obtained  by  adding  higher  order  derivatives.  The 
resulting  Taylor  Series  expansion  has  the  form 

f(U)  |u0  +  AU  "  f(Uo >  +  9f  (Up)  AU  +  1_  82f  (Uo}  (AU)2  +  1_  g3f  (Uo>  (AU)3 

3U  21  3U2  31  au3 

+  •  •  •  +  1  3  f  ^Uo  ^  (AU)n  (4.86) 

nl  au" 

If  we  make  AU  smaller,  our  accuracy  will  increase  and  higher  order  terms  can 
be  neglected  without  significant  error.  Also  since  AU  is  small,  (AU)2,  (AU)3, 
(AU)  are  very  small.  Therefore,  for  small  perturbed  values  of  U,  the 
function  can  be  accurately  approximated  by 


+  AU 


(4.87) 
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We  can  now  linearize  the  aerodynamic  forces  and  moments  using  this  tech¬ 
nique.  To  illustrate,  recall  the  lift  term  from  the  longitudinal  set  of  equa¬ 
tions  (L  =  — Fz  ) .  From  equation  4.84  we  saw  that  lift  was  a  function  of  U,  a, 
a,  Q,  .  The  Taylor  Series  expansion  for  lift  is  therefore 


L  = 


3L 


1  32L 


+  —  AU  + - AU2  + - 

3U  2  3U2 


ar  1  32h 

+  —  Aa  + - Aa2  + - 

3a  2  3a2 


+  —m  Aa  + - 

3a 


+  |Aq  + - 


+  1 3  A5e  + - 


(4.88) 


where  LQ  =  L  (U0  ,  aD  ,  ^  ,  Q0 ,  Se  ) 

o 


In  small  perturbation  theory,  each  of  the  variables  is  expressed  as  the  sum  of 
an  initial  value  plus  a  small  perturbated  value.  For  example 


and 


Therefore 


and 


U  =  U0  +  u,  where  u  =  AU  =  U  -  UQ 

0 

9u  =  9<U-V  = 

3U  3U 


3L  =  3L  3u  =  3L 
3U  3u  3U  3u 


AU  =  u 


The  second  term  of  the  expression  in  equation  4.88  then  becomes 


(4.89) 


(4.90) 


(4.91) 
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(4.92) 


9L 

au 


AU 


9L 

9u 


u 


Similarly 


9L 

9Q  AQ 


3L 

9q  q 


(4.93) 


And  all  other  terms  follow.  We  also  elect  to  let  a  =  Aa,  a  =  Aa  and  8  =  A5  . 

e  e 

Dropping  higher  order  terms  involving  u2 ,  q2 ,  etc..  Equation  4.88  now  becomes 


r  r  „  .  9L  ^  9L  •  9L  „  J  9L  t 

L  =  L_  +  —  u  +  —  a  +  —  a  +  -5-  q  +  -^  8 

9i  3<3  9Se  6 


9u 


9a 


(4.94) 


Lateral-directional  motion  is  a  function  of  8,  (3,  p,  R,  6a ,  Sr  and  can  be 
handled  in  a  similar  manner.  For  example ,  the  aerodynamic  terms  for  rolling 
moment  become 


9L  9L  9L  9La  9L, 

- 0  i - 0  4.  - p  4-  - r  +  - 6a 

98  dh  9p  9r  98a 


(4.95) 


This  development  can  be  applied  to  all  of  the  aerodynamic  forces  and  moments. 
The  equations  are  linear  and  account  for  all  variables  that  have  a  significant 
effect  on  the  aerodynamic  forces  and  moments  on  an  aircraft. 

The  equations  resulting  from  this  development  can  now  be  substituted  into 
the  LHS  of  the  equations  of  motion. 


4.6.5  Direct  Thrust  Forces  and  Moments 

Since  thrust  does  not  always  pass  through  the  eg,  its  effects  on  both  the 
force  and  moment  equations  must  be  considered  (Figure  4.25). 

The  component  of  the  thrust  vector  along  the  x-axis  is 


X^,  =  T  cos  s 

The  component  of  the  thrust  vector  along  the  z-axis  is 

Z  =  -T  sin  s 

T 


(4.96) 


(4.97) 
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FIGURE  4.25.  ORIGIN  OF  WEIGHT  AND  THRUST  EFFECTS 
ON  FORCES  AND  MOMENTS 

The  pitching  moment  component  is, 

Mr  -  T  (Zk )  =  T  Zk  (4.98) 

where  Zk  is  the  perpendicular  distance  from  the  thrust  line  to  the  eg  and  e  is 
the  thrust  angle.  For  small  disturbances,  changes  in  thrust  depend  only  upon 
the  change  in  forward  speed  and  engine  RPM.  Therefore,  by  the  same  small 
perturbation  analysis  used  for  the  aerodynamic  forces 


T  = 

T  <u'  w 

(4.99) 

T  = 

T  4.  8T  „  J. 

0  9u  u  + 

ST  X 

36  6rpk 

RPK 

(4.100) 

Thrust  effects  will  be  considered  in  the  longitudinal  equations  only  since  the 
thrust  vector  is  normally  in  the  vertical  plane  of  symmetry  and  does  not 
affect  the  lateral-directional  motion.  When  considering  engine-out 
characteristics  in  multi-engine  aircraft,  however,  the  asymmetric  thrust 
effects  must  be  considered.  Once  again,  for  clarity,  \  and  ZT  will  be 
referred  to  as  "drag  due  to  thrust”  and  "lift  due  to  thrust"  («  =  0 
assumption  in  order  to  use  the  stability  axes).  They  are  components  of  thrust 
in  the  drag  (x)  and  lift  (z)  directions.  Thus: 


4.39 


(4.101) 


*t  (To  +  au  u  +  as  5rpm)  ^cos  £) 

RPM 

3T  8T 

~  ^  +  guu  +  "95  ^rpm  ^  (sin  e)  (4.102) 

RPM 

-  (To  +  Is  u  +  |^pM  Snm>  (zk  >  (4.103) 

4.6.6  Gravity  Forces 

Gravity  acts  through  the  eg  of  an  aircraft  and,  as  a  result,  has  no 
effect  on  the  aircraft  moments.  It  does  affect  the  force  equations  as  shown 
in  Figure  4.25.  For  longitudinal  motion,  the  only  variable  to  consider  is  6. 
For  example ,  consider  the  effect  of  weight  on  the  x-axis. 

Xg  =  -mg  sin  6  (4.104) 

Since  m  and  g  are  considered  constant,  0  is  the  only  pertinent  variable. 
Therefore,  the  expansion  of  the  gravity  term,  Xg ,  can  be  expressed  using  the 
small  perturbation  assumption  as 

3X 

xq  =  Xq  + _ I  6  (X  =  equilibrium  condition  of  X  )  (4.105) 

yo  ae  9o  g 

For  simplification,  the  term  Xg  will  be  referred  to  as  drag  due  to 
weight,  (Dwt)-  This  incorporates  the  small  angle  assumption  that  was  made  in 
development  of  the  aerodynamic  terms;  however,  the  effect  is  negligible. 
Therefore,  Equation  4.105  becomes 

3D 

Dwt  =  Do  +  e  (4.io6) 

wt  ae 

Likewise  the  z-force  can  be  expressed  as  negative  lift  due  to  weight  (Lwt), 
and  the  expanded  term  becomes 

3L 

Lwt  =  L0  +  6  (4.107) 

ae 

The  effect  of  gravity  on  side  force  depends  solely  on  bank  angle  ( <M  ,assuminq 
small  0.  Therefore,  y 


4.40 


(4.108) 


=  Y0  +  ♦ 

wt  9(J> 

These  component  equations  relate  the  effects  of  gravity  to  the  equations  of 
motion  and  can  be  substituted  into  the  LHS  of  the  equations. 

4.6.7  Gyroscopic  Moments 

Gyroscopic  effects  are  insignificant  for  most  static  and  dynamic  analyses 
since  angular  rates  are  not  considered  large.  They  begin  to  become  important 
as  angular  rates  increase  (i.e.,  P,  Q,  and  R  become  large).  For  spin  and  roll 
coupling  analyses,  they  are  large  and  gyroscopic  effects  must  be  considered. 
In  this  basic  development  of  the  equations  of  motion,  however,  they  will  be 
assumed  to  be  negligible. 

4.6.8  Expanded  LHS  Equations 

Using  the  previous  developments,  the  expanded  LHS  equations  become 


"DRAG"  -[D„  +|eu+|5a  +  i|i  +  |?q+  ®  5J  +  [To  +^u  +  ||_S,t„]  (CCS  e) 


RPM 


-  +  IN 


(4.109) 


wt 


■LIFT-  -IL„  ♦  ||u  +  |ia  +  Si  +  ||q  ♦  %  «.  I  ♦  [-T,  +  ]  (sin  ,) 

da  ^  e  RPM 


+  +  %e] 


(4.110) 


"YAW"  N 


9MA 

3M 

A 

8V 

3M 

A 

9Ma  ,  , 

+  [T0 

+ - u  + 

- a  + 

a 

+  - q 

+  -  «.] 

3u 

3a 

3a 

3q 

3Se 

3Ya 

SYa. 

3Y, 

3Y, 

3Y. 

9Ya 

+  w*  + 

——6  +  ■ 
33 

3pP+  9t  r+  3S,5.  + 

3T5.  + 

r 

3La 

3La. 

3La 

3La 

9La 

3L 

+  3Te  + 

— e  + 

33 

3P"P 

+  3r~r 

+  3S~Sa  + 

a 

*T«r 

r 

3N 

3N 

3Na 

3Na 

3N 

3N 

+  3TP  +  - 

— 3  + 
33 

3p"P 

+  3r“r 

+  inr*  + 
36  a 

a 

36  *r 

r 

3T  .  3T 


RPM 


3Y 

\,+  * 


(4.111) 


(4.113) 


:4.114) 
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4.7  RHS  IN  TERMS  OF  SMALL  PERTURBATIONS 


To  conform  with  the  Taylor  Series  expansion  of  the  LHS,  the  RHS  must  also 
be  expressed  in  terms  of  small  perturbations.  Recall  that  each  variable  is 

expressed  as  the  sum  of  an  equilibrium  value  plus  a  small  perturbed  value 

(i.e.,  U  =  U0  +  u,  Q  =  Q0  +  q,  etc.).  These  expressions  can  be  substituted 
directly  into  the  full  set  of  the  RHS  equations  (Equations  4.23  -  4.25,  and 
4.59  -  4.61).  As  an  example,  the  lift  equation  (z-direction  of  longitudinal 
equations)  will  be  expanded.  Start  with  the  RHS  of  Equation  4.25 

Fj.  -  m  (W  +  PV  -  QU)  (4.25) 

Substitute  the  initial  plus  perturbed  values  for  each  variable. 

Fz  =  m  [W0  +  w  +  (P0  +  p)  (V0  +  v)  -  (Q0  +  q)  (Up  +  u)]  (4.116) 

Multiplying  out  each  term  yields: 

Fz  =  m  [W0  +  w  +  P0  V0  +  p  V0  +  P0v  +  pv  -  Q0  Up  -  q  Up  -  Qqu  -  qu]  (4.117) 

Applying  the  boundary  conditions,  (assumptions  from  subsection  4. 6. 4.3. 3), 

simplifies  the  equation  to 


Fz  =  m  [w  +  pv  -  qU0  -  qu]  «  m  [w  +  pv  -  q  (Un  +  u)]  (4.118) 

or 


F 

Z 

=  m 

[w 

+  pv  -  qU] 

(4, 

.119) 

Using  this  same  technique. 

the  ! 

set 

of  RHS  equations  become: 

Longitudinal 

"DRAG" : 

m 

(u  + 

qw 

-  rv) 

(4. 

.120) 

"LIFT" : 

m 

(w  + 

pv 

-  qU) 

(4. 

.121) 

"PITCH" : 

• 

q 

I  - 

y 

Pr 

U.  -  I.)  +  (P2  -  r2)  i„ 

(4, 

,122) 

Lateral-Directional 


'SIDE": 

m  (v  +  rU 

-  pw) 

(4.123) 

ROLL" : 

P  Ix  +  qr 

-  y 

-  (r  +  pq)  Ixz 

(4.124) 

YAW": 

r  \  +  pq 

(i  - 1 ) 

y  x 

+  (qr  -  p  )  I„ 

(4.125) 

Which  are  valid  for  small  perturbations  about  an  equilibrium  condition  of 
steady  straight  symmetric  flight. 


4.8  REDUCTION  OF  EQUATIONS  TO  A  USABLE  FORM 


4.8.1  Normalization  Of  Equations 

To  put  the  linearized  expressions  into  a  more  usable  form,  each  equation 
is  multiplied  by  a  "normalization  factor."  This  factor  is  different  for  each 
equation  and  is  picked  to  simplify  the  first  term  on  the  RHS  of  the  equation. 
It  is  desirable  to  have  the  first  term  of  the  RHS  be  either  a  pure 
acceleration  (u,  p,  q,  or  r),  or  angular  rate  (a,  or  (3)  and  these  terms  were 
previously  identified  in  equations  4.84  and  4.85  as  the  longitudinal  or 
lateral-directional  variables.  As  shown  previously 


and 


Since  we  have  assumed  that  V„ 

T 


and 


-  U0  and  V0 

•  w 


«  W„ 


=  0,  we  have 


(4.16) 

(4.18) 


(4.126) 

(4.127) 


Table  4.1  shows  the  normalizing  factors  and  subsequent  equations  of  motion. 


TABLE  4.1 


NORMALIZING  FACTORS 


Equation 


Normalizing 

Factor 


First  Term  is  Now 
Pure  Accel/Ang  Rate 


Units 


"DRAG" 

1 

-2  + 

1 

+ 

*rl 

• 

-  =  u  - 

— 

t  £t2  1 

(4.128 

m 

m 

m 

sec 

"LIFT" 

1 

mu0 

L 

mUQ 

2t 

+  iu0+- 

- - = 

a  +  - 

■ 

secJ 
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4.8.2  Stability  Parameters 

Stability  parameters  are  quantities  that  express  the  variation  of  force 
or  moment  on  the  aircraft  caused  by  a  distrurbance  from  steady  flight.  They 
are  simply  the  partial  coefficients  (9L/9u,  etc.)  multiplied  by  their 
respective  normalizing  factors.  They  express  the  variation  of  forces  or 
moments  caused  by  a  disturbance  from  steady  state.  Stability  parameters  are 
important  because  they  can  be  used  directly  as  numerical  coefficients  in  a  set 
of  simultaneous  differential  equations  describing  the  dynamics  of  an  airframe. 
To  demonstrate  their  development,  consider  the  aerodynamic  terms  of  the  lift 
equation.  By  multiplying  Equation  4.94  by  the  normalizing  factor  l/mU0 ,  we  get 


+  k  +  A £a  +  A  fii  +  +  A **  [li](4-134) 
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The  indicated  quantities  are  defined  as  stability  parameters  and  the 
equation  becomes 


jrad] 

LsecJ 


(4.135) 


Stability  parameters  have  various  dimensions  depending  on  whether  they  are 
multiplied  by  a  linear  velocity,  an  angle,  or  an  angular  rate.  For  example, 
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The  lateral-directional  motion  can  be  handled  in  a  similar  manner.  For 
example,  the  normalized  aerodynamic  rolling  moment  becomes 


L  L  L  L  •  L  L  L  *  L  * 
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sec 


These  stability  parameters  are  sometimes  called  "dimensional  derivatives"  or 
"stability  derivative  parameters,"  but  we  will  reserve  the  word  "derivative" 
to  indicate  the  nondimensional  form  which  can  be  obtained  by  rearrangement. 
This  will  be  developed  later  in  the  chapter.  See  subsection  4.9  for  a 
complete  set  of  equations  in  stability  parameter  form. 


4.8.3  Simplification  Of  The  Equations 

By  combining  all  of  the  terms  derived  so  far,  the  resulting  equations  are 
somewhat  lengthy.  In  order  to  economize  on  effort,  several  simplifications 
can  be  made.  For  one,  all  "small  effect"  terms  can  be  disregarded.  Normally 
these  terms  are  an  order  of  magnitude  less  than  the  more  predominant  terms. 
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These  and  other  simplifications  will  help  derive  a  concise  and  workable  set  of 
equations. 


4.8.4  Longitudinal  Equations 


4. 8. 4.1  Drag  Equation.  The  complete  normalized  drag  equation  is 


Aero  Terms 
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Simplifying  assumptions 


T  D  Do 

1.  o  cos  e  -  o  -  wt  ~  o 


m 


m  m 


(Steady  State,  Sum  to  Zero) 


2  ri?  „  +  3T 
’  l3u  U  38, 


5epm](co5  e)  =  0  (Constant  RPM,  is  small) 


RPM 


3u 


3.  rv  =  0 


(No  lat-dir  motion) 

The  small  perturbation  assumption 
allows  us  to  analyze  the  longitudinal 
motion  independent  of  lateral- 
directional  motion. 


4.  qw  si  0 


(Order  of  magnitude) 


5.  D^f  and  Da  are  all  very  small,  essentially  zero. 


The  resulting  equation  is 


-  [Daa  +  Duu  +  De0  +  D5  Se  ]  =  u 


(4.138) 
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Rearranging 


-  [u  +  Duu  +  Daa+Dee]  =  D&  Se 

e 

4.8. 4.2  Lift  Equation.  The  complete  lift  equation  is 


(4.139) 
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Simplifying  assumptions 
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2. 
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(Constant  RPM, -5-  is  small) 

dU 

(Order  of  magnitude,  for 
small  0) 


4.  pv  =  0 


(No  lat-dir  motion) 


5. 


sy  «  q 
u. 


(U  =  U0) 


The  resulting  equation  is 

-  [L  a  +  L.a  +  Lu  +  Lq  +  L.  6]  = 

a  a  u  q  ~  6  e  J 
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Rearranging 


a  -  q 


(4.141) 
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4. 8. 4. 3  Pitch  Moment  Equation. 


_ 1  +  Ma  “  +  +  MA  u  +  Ma  $e  +  Ma  q  +  thrust  Terms 

Iy  u  6  ^  q 

=  q  -  pr  [ ^  +  (P  “  r  ^  i 


(4.143) 


This  can  be  simplified  as  before  (using  the  simplifying  assumptions).  NOTE: 
For  convenience,  we  will  also  drop  the  "A"  subscripts.  Thus 

q  -  Mqq  -  Maoc  -  M&a  -  Muu  =  M  &e  (4.144) 

& 

e 

We  now  have  three  longitudinal  equations  that  are  easy  to  work  with.  Notice 
that  there  are  four  variables,  9,  a,  u,  and  q,  but  only  three  equations.  To 
solve  this  problem,  0  can  be  substituted  for  q. 


q  =  9,  and  q  =  0 


This  can  be  verified  from  the  Euler  angle  transformation  for  pitch  rate  where 
the  roll  angle,  <f>,  is  zero. 


q 


<f>  +  ^  zy? < 


0 

<(>  cos  9 


.'.  q  =  0 

This  gives  us  another  restriction  to  this  development  of  the  equations  of 
motion: 


ASSUMPTION: 


The  initial  bank  angle  is  zero,  <f>0  =  0. 
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4.8.5  Lateral-Directional  Equations 

The  complete  lateral-directional  equations  are  as  follows: 


J±  +  \e*\e*\p*\*  +  V  «.  +  ^ 

mlL  P  •  P  r  &  & 

u  P  a  r 

+  °wt+Yi<f>  = 
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(4.145) 


+  La  P  +  La.8  +  La  p  +  La  r  +  La  5a  +  La  6r 

P  8  P  r  6  6 

a  r 


=  p  +  qr 


(r  +  pq) 


(4.146) 


\  +  NA  3  +  +  NA  P  +  NA  r  +  NA  Sa  +  NA  Sr 

r~  p  e  P  r  6a  6r 

"  r  +  pq  ( Iy  ~  )  +  (qr  -  p)  (4.147) 


In  order  to  simplify  the  equations ,  the  following  assumptions  are  made: 

1.  A  wings  level  steady  state  condition  exists  initially.  Therefore, 

La  ,  Na  ,  Ya  ,  and  Y0  are  zero. 

0  0  0  wt 

2.  p  =  <f>,  p  =  4>  (0  =  0,  see  Euler  angle  transforma¬ 

tions  for  roll  rate,  equation  4.70) 


3.  The  terms  YA  8,  LA  8/  and  NA  8  are  all  small,  essentially  zero. 
p  p  p 


4. 


5.  q  s  0 

6.  pw  =  0 


(no  pitching  rate  in  the  Lateral- 
Direction  equations) 

(no  longitudinal  motion) 


NOTE:  We  will  also  drop  the  "A"  subscripts  for  convenience. 


4.9  STABILITY  PARAMETER  FORM  OF  THE  EQUATIONS  OF  MOTION 

Using  the  previous  subsection,  the  equations  of  motion  reduce  to  the 
following  stability  parameter  form: 
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Lateral 

Directional  Equations 

(3) 

(<f>) 

(r) 

"SIDE 

FORCE" 

e-V 

- 

+  (1-Yr)  r 

ya  s  +  y..  s 

6  a  6  r 

a  r 

(4.151) 
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"YAWING" 
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There  are  six  unknowns  and  six  equations.  The  terms  on  the  RHS  are  now  the 
inputs  or  "forcing  functions."  Therefore,  for  any  input  Se ,  8a ,  or  S  ,  the 
equations  can  be  used  solved  to  get  0,  u,  a,  0,  <j>,  and  r. 

Stability  parameter  data  is  usually  the  format  found  in  textbooks  or  NASA 
published  aircraft  stability  and  control  data  reports. 

4.10  STABILITY  DERIVATIVES 

The  parametric  equations  give  all  the  information  necessary  to  describe 
the  motion  of  any  particular  aircraft.  There  is  only  one  problem.  When  using 
a  wind  tunnel  model  for  verification,  a  scaling  factor  must  be  used  to  find 
the  values  for  the  aircraft.  It  is  difficult,  therefore,  to  compare  aircraft 
using  stability  parameters.  In  order  to  eliminate  this  requirement,  a  set  of 
nondimensional  equations  can  be  derived.  This  is  best  illustrated  by  an 
example : 

Given  the  parametric  equation  for  pitching  moment, 

e  _  Mqe  -  Muu  -  Maa  -  =  Ms  Se  (4.154) 

e 

Let  us  derive  an  equation  in  which  all  terms  are  NONDIMENSIONAL . 

The  steps  in  this  process  are: 

1.  Take  each  stability  parameter  and  substitute  its  coefficient  relation 
and  take  the  derivative  at  the  initial  condition,  with  other 
variables  (i.e.  u,  a,  a,  5e )  held  constant.  From  subsonic 
aerodynamics, 

M  -  Cm  -  pl^Sc  =  Mqq  +  MuU  +  M.i  +  M^a  -  Se  (4.155) 

2  e 

Looking  at  the  change  in  this  pitching  moment  due  to  a  change  in  pitch  rate 
(as  an  example)  at  the  initial  conditions 

3(C  1  p  U2 Sc) 

Ma  =  i_  2 _  (4.156) 

!y  aq  iy  aq  o 
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(4.157) 


Cm  is  the  only  variable  that  is  dependent  on  q,  therefore, 

H  =  pUo2  Sc  9Cm 
9  2ly  3q 


2.  Nondimensionalize  the  partial  term: 

3C 

w  ^ “-**»  - 

To  nondimensionalize  the  partial  terms,  certain  compensating  factors  are 
customarily  used  (Table  4.2).  in  this  case,  the  compensating  factor  is 


_c_  [ft] 
2U0  [ft/sec] 


sec 


Multiply  and  divide  Equation  (4.157)  by  the  compensating  factor  and  get 
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<3 


Checking 


PU0  Sc  2U0 


21 


ac 
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aq 


>  =>  This  term  is  dimensionless  (4.158) 


dimensionless 

ft/sec 

ft/sec 


dimensionless 


This  is  called  a  stability  derivative  and  is  written 


ac 


The  basic  nondimensional  form  Cm  is  important  because  correlation  be- 

q 

tween  geometrically  similar  airframes  or  the  same  airframe  at  different  flight 
conditions  is  easily  attained  with  stability  derivatives  (this  cannot  easily 
be  done  with  stability  parameters).  Additionally,  aerodynamic  stability 
derivative  data  from  wind  tunnel  tests,  flight  tests,  and  theoretical  analyses 
are  usually  presented  in  nondimensional  form. 

Stability  derivatives  generally  fall  into  two  classes:  static  and 

dynamic.  Static  derivatives  arise  from  the  position  of  the  airframe  with 
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respect  to  the  relative  wind  (i.e.,  CL  ,  C,  ,  Cn  ,  q  ).  Whereas,  dynamic 
derivatives  arise  from  the  motion  (velocities)* of  tfie  airframe  i.e.,  C  ,  C  , 

cn  ,  cx  ).  “  u 

r  P 

3.  The  entire  term  with  stability  derivative  is 

pu02  Sc 

Mqq  21  2U  q  (4.160) 

y  0  q 

4.  We  can  do  the  same  for  each  term  in  the  parametric  equation.  For 

example , 
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Since  both  Cm  and  U  are  functions  of  u,  then 
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0  since  the  initial  conditions  are  steady  state.  The  compensating 
factor  for  this  case  is  1/UQ 


Muu 


PU0'  Sc 
21 


u 
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5.  Once  all  of  the  terms  have  been  derived,  they  are  substituted  into 
the  original  equation,  and  multiplied  by 
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compensating  factors 

for  all  of  the  variables  are 

listed  in  Table 

4.2. 

Variable 

TABLE  4.2 

COMPENSATING  FACTORS 
Compensating 

Factor 

Nondimensi onal 
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II 
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a 

(rad/sec) 
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2U0 

Cot 

^0 
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(ft/sec) 

1 

Uo 

u 

Uo 

a 

( radians ) 
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a 

e 

( radians) 
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e 

4.10  STABILITY  DERIVATIVE  FORM  OF  THE  EQUATIONS  OF  MOTION 


The  simplified  equations  of  motion  in  stability  derivative  form  are  shown 
in  Table  4.3.  The  derivation  of  these  equations  has  been  presented  to  give  an 
understanding  of  their  origin  and  what  they  represent.  It  is  not  necessary  to 
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LONGITUDINAL  EQUATIONS 


be  able  to  derive  each  and  every  one  of  the  equations.  It  is  important, 
however,  to  understand  several  facts  about  the  nondimensional  equations. 


1.  Since  these  equations  are  nondimensional,  they  can  be  used  to 
compare  aircraft  characteristics  of  geometrically  similar  airframes. 

2.  Stability  derivatives  can  be  thought  of  as  if  they  were  stability 
parameters.  Therefore,  Cm  refers  to  the  same  aerodynamic 

a 

characteristics  as  Ma ,  only  it  is  in  a  nondimensional  form. 

3.  Most  aircraft  designers  and  builders  are  accustomed  to  speaking  in 
terms  of  stability  derivatives.  Therefore,  it  is  a  good  idea  to 
develop  a  "feel"  for  all  of  the  important  ones. 

4.  These  equations  as  well  as  the  parametric  equations  (eq.  4.70  to 

4.72)  describe  the  complete  motion  of  an  aircraft.  They  can  be 

programmed  directly  into  a  computer  and  connected  to  a  flight 
simulator.  They  may  also  be  used  in  cursory  design  analyses.  Due 
to  their  simplicity,  they  are  especially  useful  as  an  analytical 
tool  to  investigate  aircraft  handling  qualities  and  determine  the 
effect  of  changes  in  aircraft  design. 

5.  Sometimes  non-dimensional  stability  derivative  data  are  tabulated  as 
a  function  of  a  and  Mach  number  to  cover  the  entire  aircraft  flight 
envelope.  Equations  of  motion  for  other  than  1  g  flight  can  be 
obtained  from  this  format  given  a  corresponding  to  a  particular  load 
factor . 

4.11  AIRCRAFT  TRANSFER  FUNCTIONS 

If  the  Laplace  transform  of  the  stability  parameter  equations  of  motion 
(subsection  4.9)  are  taken,  assuming  zero  initial  conditions,  and  the 
equations  are  written  in  matrix  notation,  the  following  equations  result  (in 
terms  of  stability  parameters  for  the  aircraft  stability  axis  system): 
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The  above  matrix  equations  are  all  that  are  needed  for  transfer  function 
derivation.  As  an  example,  the  transfer  function  relating  a  to  S  is 


(4.176) 


which  can  be  found  from  Cramer's  rule  as 
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where  the  determinant  in  the  denominator  yields  the  characteristic  equation  of 
the  unaugmented  aircraft.  Similar  derivations  can  be  used  to  find  any  single 
input,  single  output  transfer  function  for  the  longitudinal  or 
lateral-directional  axes. 
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4 -12  AIRCRAFT  TRANSFER  FUNCTION  DATA  PRESENTATION 


Examples  of  three  formats  used  to  present  aircraft  transfer  function  data 
are  shown  in  Figure  4.26  and  Tables  4.4  and  4.5  for  the  lateral-directional 
axes  of  the  A-7A. 

Hie  data  presented  in  Figure  4.26  can  be  used  in  the  non-dimensional 
stability  derivative  equations.  The  data  presented  in  Table  4.4  can  be  used 
in  the  dimensional  stability  parameter  equations.  Both  approaches  yield  the 
same  transfer  functions.  The  data  in  Table  4.5  are  presented  in  transfer 
function  format  already  (body  axis  system  in  this  case).  The  symbol  A  is  the 
lateral-directional  characteristic  equation  (denominator  term).  The  numerator 
terms  are  presented  as 


..output 

input 


(s) 


for  example  (S) 

a 

is  the  numerator  term  for  roll  rate  due  to  aileron  deflection.  1/T  denotes 
real  axis  poles  or  zeros  (discussed  in  Chapter  13).  C  and  w  denotes  damping 
ratio  and  natural  frequency  of  poles  or  zeros  (regardless  of  the  subscript  in 
the  table).  These  must  be  converted  to  the  form 


s  =  -  <r  +  wdj  (4.178) 

where 

a  =  tcon  and  &>d  =  wn  /  1  -  C2  (4.179) 

A  is  the  root  locus  gain  (discussed  in  Chapter  13)  of  the  transfer  function. 
For  example,  the  transfer  function 


G*  (s)  = 


K  (s) 


A  (s) 


(4.180) 


at  0.6  Mach,  sea  level  has  the  characteristic  equation 


Ms)  =  (s  +  0.0411)  (s  +  4.46)  [s2  +  2(0.202)  (2. 91)s  +  (2.912)]  (4.181) 

or 

Ms)  m  (s  +  0.0411)  ( s  +  4.46)  [s  +  0.59  ±  2.85j)  (4.182) 

and 

(s)  =  28.4 (s  -  0.00234) (s2  +  2(0.217) (3. 05)s  +  (3.05)2)  (4.183) 

a 

The  units  are  radians  for  &a  and  radians  per  second  for  p,  or  degrees  for  Sa 
and  degrees  per  second  for  p,  since  both  the  input  and  output  parameters  have 
similar  units.  If  the  input  and  output  parameters  have  different  units,  as  is 
the  case  in  the  transfer  function 

N/(s) 

a 

then  the  units  are  feet  per  second  squared  for  ay  and  radians  for  5a .  If  in 
doubt,  always  assume  radians  for  angular  parameters. 

Tables  4.6  and  4.7  give  typical  values  for  stability  derivatives.  Note 
that  the  tables  use  the  general  body  axes  (not  the  stability  axes)  which  use 
the  forces  X  and  Z,  not  L  and  D. 


0 


.2 


.4 


FIGURE  4.26  TYPICAL  STABILITY  DERIVATIVE  DATA  (continued) 
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6,=  -^5 - (INCLUDES  SPOILER  EFFECTS) 


FIGURE  4.26  TYPICAL  STABILITY  DERIVATIVE  DATA  (continued) 


4.63 


FIGURE  4.26  TYPICAL  STABILITY  DERIVATIVE  DATA  (continued) 
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TABLE  4.4 

LATERAL-DIRECTIONAL  STABILITY  PARAMETERS  FOR  THE  A-7A 


Note:  Data  are  for  body-fixed  centerline  axes,  dean  flexible  airplane 


■ 

FUGHT  CONDITION 

■ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

0 

0 

0 

15000 

15000 

15O00 

15000 

35000 

35000 

M 

0.25 

0.6 

03 

03 

03 

03 

LI 

03 

03 

Y* 

-0.162 

-0314 

-0314 

-0.122 

-0.187 

-0310 

—0435 

-00847 

-0.145 

Y6. 

-0.00274 

-0.0105 

-0.00857 

-000150 

-000655 

-0.00691 

-000216 

-000267 

— 0j004J7 

Y5, 

0.0430 

0.0769 

0.0626 

00307 

0.0537 

00550 

00192 

00267 

00347 

10* 

-119 

-443 

—98.0 

-8.79 

-293 

-660 

-713 

-143 

-303 

H 

-2.00 

—4.46 

-9.75 

-138 

-2.73 

-6.19 

-731 

-140 

-300 

1.18 

1.15 

138 

0357 

0368 

0343 

0869 

0399 

0363 

L6. 

534 

28.4 

253 

3.75* 

173 

24.1 

123 

736 

143 

232 

1L4 

133 

132 

737 

113 

737 

309 

635 

n3# 

138 

5.74 

173 

0348 

£12 

103 

213 

138 

4.72 

NP 

“0.0870 

-0.168 

-0319 

-0.0310 

-0.116 

-0307 

-0.169 

-0.0799 

-0.112 

Nr 

-0369 

-0.905 

-134 

-0371 

-0341 

-0375 

-133 

-0347 

-0455 

N6. 

0.402 

208 

136 

0380 

137 

1.64 

L04 

0352 

101 

% 

-133 

-8.61 

— 1L1 

-136 

-534 

-830 

-433 

-234 

-5.11 

*S  *  v/vT 


TABLE  4.5 

AILERON  LATERAL  TRANSFER  FUNCTION  FACTORS  FOR  THE  A-7A 


Note:  Data  are  for  body-fixed  centerline  axes,  dean  flexible  airplane 


FUGHT  CONDITION 

1 

2 

3 

4 

5 

6 

7 

8 

9 

h 

0 

n 

0 

15.000 

15,000 

15,000 

15,000 

35,000 

M 

0.25 

03 

03 

0.6 

03 

LI 

0.6 

03 

m. 

0.0462 

0.0411 

0.0180 

0.0449 

0.0435 

0.0214 

0.0102 

0.0319 

0.0191 

A 

1/Tr 

L62 

4.46 

9.75 

0368 

2.71 

6.17 

7.15 

128 

292 

?d 

0237 

0202 

0218 

0231 

0.156 

0.175 

0.189 

0.114 

0.128 

«d 

131 

231 

4.68 

1.65 

229 

3.66 

5.03 

131 

258 

El 

Ap 

534 

28.4 

252 

325 

17.6 

24.1 

122 

736 

142 

H 

—0.0219 

-0.00234 

-0.00113 

-0.0232 

-0.00347 

-0.00144 

-0.00137 

-0.00718 

-0.00241 

0217 

0217 

0222 

0.191 

0.173 

0.176 

0.173 

0.122 

0.124 

■ 

«p 

1.49 

aos 

4.91 

127 

224 

337 

523 

L62 

264 

<t> 

N«. 

A</> 

5.42 

283 

252 

331 

17.7 

24.1 

125 

8.04 

142 

4 

0210 

0217 

0222 

0.183 

0.173 

0.177 

0.175 

0.119 

0.124 

a<t> 

LSI 

aos 

431 

129 

224 

337 

522 

1.62 

264 

N6 

A, 

0.402 

2.08 

136 

0280 

127 

1.64 

1.04 

0.652 

L01 

0396 

L12 

1.13 

0.445 

0.777 

0.944 

0281 

0.420 

0293 

f. 

0.0852 

0287 

0397 

0.146 

0.151 

0.446 

0.638 

0.0198 

0.193 

“>r 

235 

229 

326 

2.18 

213 

278 

298 

203 

245 

4. 

A0 

-0.00274 

-0.0105 

-0.00857 

-0.00150 

-0.00655 

-0.00691 

-0.00216 

-0.00267 

-0.00427 

(0.885) 

326 

7.76 

(0.726) 

221 

5.77 

10.7 

0.793 

(0372) 

(0.667) 

-0.627 

-0254 

(0.471) 

-133 

-0245 

-0.113 

-0.422 

(102) 

-233 

63.1 

782 

-391 

232 

863 

188 

-147 

-0245 

< 

\ 

-0.766 

-7.06 

-8.61 

-0.477 

-4.16 

—628 

-221 

-126 

-0.0374 

CG 

"VS? 

(0.943) 

229 

-1.16 

(0.758) 

122 

-0296 

-0.146 

0290 

(0.801) 

1/T.  («,  > 

'2  S 

(0.648) 

532 

-134 

(0.461) 

3.12 

-266 

-733 

0361 

(234) 

V1^ 

0.0696 

—0310 

(3.65) 

0.0673 

-0294 

(279) 

0397 

0.0499 

-0.113 

637 

L76 

(10.7) 

7.10 

139 

(—633) 

931 

332 

120 

TABLE  4.5  (concluded) 

AILERON  LATERAL  TRANSFER  FUNCTION  FACTORS  FOR  TOE  A- 
Note:  Data  are  for  body-fixed  centerline  axeidean  flexible  airplane 


FLIGHT  CONDITION 


»«,  »6 


-l  \ 

CG  l/T 


-0.0624  -0.00199  0.00026 6  -0.0603  -0.00616  a000578  a00271  I  -00178  -a00216 


-0.0145  -0.00502  -0.0227  -0.00654 


-00436  I  -0.0107 


TABLE  4.6 

TYPICAL  VALUES  FOR  LONGITUDINAL  STABILITY  DERIVATIVES 


STABILITY 

DERIVATIVE 

EQUATION 

TYPICAL 

VALUE 

CV 

“2Cd  -  U0  l  3Cd/  flu] 

-0.05 

C*« 

CL  -  [  dCD  /  d  <*  J 

+0.1 

C*u 

“2Cl  -  U0  [  SCL/  du] 

-.05 

C*U 

— —  [  dCL/  ft  «  ] 

-4 

C*ci 

-1 

s 

-2 

Cm„ 

Neglect  for  Jets 

Cma 

-0.3 

C"a 

-3 

C"q 

-8 

NOTE:  EFFECTS  OF  Cm^ANDCm<|  ARE  USUALLY  COMBINED  WHEN  USING 
FLIGHT  TEST  DATA. 


TABLE  4.7 

TYPICAL  VALUES  FOR  LATERAL-DIRECTIONAL  STABILITY  DERIVATIVES 

STABILITY  TYPICAL 

DERIVATIVE  VALUE 


s 

-0.6 

s 

-0.06 

c'p 

-0.4 

c', 

+0.06 

% 

+0.11 

c"p 

-0.015 

C"r 

-0.12 

s 

Neglect 

Neglect 
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4.13  LIST  OF  ABBREVIATIONS  AND  SYMBOLS 


Symbol  Definition 

A  Gain  of  root  locus  transfer  function 


b 


Wingspan 


c  Mean  Aerodynamic  Chord:  The  theoretical  chord  for  a  wing  which 

has  the  same  force  vector  as  the  actual  wing  (also  MAC). 

c9  Center  of  Gravity:  common  name  for  the  aircraft's  center  of 

mass. 


D 


Wt 


wt 


f 


F 


Drag:  The  component  of  the  resultant  aerodynamic  force 

parallel  to  the  relative  wind.  It  must  be  specified 
whether  this  applies  to  a  complete  aircraft  or  to  parts 
thereof. 

Drag  due  to  weight. 

Initial  drag  due  to  weight. 

function 

Applied  force  vector. 


g 


G  ,G  , G 

x  9  y  f  z 


Components  of  applied  forces  on  respective  body  axes. 

Local  gravitational  vector  where  g  is  the  magnitude  of  the 
acceleration  due  to  gravity. 

Applied  moment  vector. 

Components  of  the  applied  moments  on  the  respective  body  axes. 


H 


Hx  'Hy  'Hz 

I,  5,  K 


I  ,1  ,1 

x  9  y  9  z 


■'"x  y  '  z  '  *x  z 


Angular  momentum  vector. 

Components  of  the  angular  momentum  vector  on  the  body  axes. 

Unit  vectors  in  the  body  axis  system. 

Moments  of  inertia  about  respective  body  axes. 

(Example:  Ix  =  Jv  pA  (y2  +  z2  )  dV) 

Products  of  inertia,  a  measure  of  symmetry. 

(Example:  Ixy  =  I  yx  =  Jv  xypAdV) 
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Symbol 


Definition 


Lift:  The  component  of  the  resultant  aerodynamic  force 
perpendicular  to  the  relative  wind.  It  must  be  specified 
whether  this  applies  to  a  complete  aircraft  or  to  parts 
thereof. 


La'  V  Na 
Loth '^oth 'N0th 

V  V 
L,  or  1 

Lwt 

Lo 

wt 

LHS 


M 


Aerodynamic  moments  about  the  x,  y  and  z  vehicle  axes. 

Other  moments  about  the  x,  y  and  z  vehicle  axes. 

Thrust  moments  about  the  x,  y  and  z  vehicle  axes. 

Rolling  moment. 

Lift  due  to  weight 
Initial  lift  due  to  weight 

Left  hand  side  -  the  side  of  the  equations  of  motion  which 
represents  the  applied  forces  and  moments  on  the  aircraft. 

Mach  number 


M,  or  m  Pitching  moment 

m  mass  of  aircraft 


N,  or  n  Yawing  Moment 


P,Q,R 


V  Qo •  Ro 

P,q/r 

r 

Re 

RPM 

RW 


Angular  rates  about  the  x,  y,  and  z  vehicle  axes. 

Initial  angular  rates. 

Perturbed  values  of  P,Q,R,  respectively. 

Position  vector  measured  from  the  eg. 

Reynolds  number. 

Revolutions  per  minute,  used  to  indicate  engine  throttle 
setting. 

Relative  wind 


Right  hand  side  -  the  side  of  the  equations  of  motion  which 
represents  the  aircraft's  response  to  the  applied  forces  and 
moments . 


s 

S 

t 


Laplace  transform  variable. 
Wing  area. 

Time 
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Symbol 

T 


To 

u,vfw 


Uo  'Vo  'Wo 


u,v,w 

X 

V 

V 

VT 

wt 

X,  Y,Z 
X,Y,Z 

WZA 

VVZ« 

X  ,Y  ,Z 

g  r  g  9  g 
0  0  0 

X  .  ,Y  .  .  ,Z 

othf  oth  f  oth 

VYt'Zt 

x,y,z 


Xs  'Ys  ,Zs 


wt 


wt 


a 


Definition 

Thrust;  temperature;  time  constant 
Initial  thrust. 

Components  of  velocity  along  the  x,  y,  and  z  vehicle  axes. 

Components  of  velocity  along  the  x,  y,  and  z  vehicle  axes  at 
zero  time  (i.e.,  equilibrium  condition). 

Small  perturbations  of  U,V,W,  respectively. 

Vector  cross  product. 

Aircraft  volume 

velocity 

True  velocity  (the  velocity  of  the  relative  wind). 

Weight. 

Axes  in  the  earth  axis  system. 

Aircraft  velocities  in  the  inertial  coordinate  system. 
Aerodynamic  forces  in  the  vehicle  axes. 

Gravity  force  in  the  vehicle  axes. 

Initial  gravity  forces. 

Other  forces  in  the  vehicle  axes. 

Thrust  forces  in  the  vehicle  axes. 

Axes  in  the  body  axis  system. 

Axes  in  the  stability  axis  system. 

Side  force  due  to  weight. 

Initial  side  force  due  to  weight. 

Distance  between  the  thrust  line  and  the  eg. 

Air  density. 

Mass  density  of  aircraft. 

Angle  of  attack,  also  a  small  perturbation  in  a. 
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Symbol 

“o 

(3 

e0 

A 

K  '  Se  '  ' 
^RPM 

C 

a 

C 

\p,e,<i> 

00 

"d 


Definition 
Initial  a. 

Sideslip  angle,  also  a  small  perturbation  in  3. 

Initial  3. 

Characteristic  equation;  also  used  to  denote  a  change  in  a 
variable. 

Deflection  angle  of  the  ailerons,  elevator,  and  rudder, 
respectively,  also  small  perturbations  in  these  values. 

Change  in  engine  throttle  setting. 

Thrust  angle. 

Total  damping,  a  =  £oon . 

Damping  ratio. 

Euler  angles:  yaw,  pitch,  and  roll,  respectively. 

Total  angular  velocity  vector  of  an  aircraft, 
damped  natural  frequency, 
natural  frequency, 
therefore. 


PROBLEMS 


4.1.  Draw  in  the  vectors  VT,  U,  V,  W,  and  show  the  angles  a  and  (3. 


Derive  a  -  W/VT  and  3  =  V/VT  using  the  small  angle  assumption. 

4.2.  Define  "right  hand"  and  "orthogonal"  with  reference  to  a  coordinate 
system. 

4.3.  Define  "moving"  earth  and  "fixed"  earth  axis  coordinate  systems. 

4.4.  Describe  the  "body"  and  "stability"  axes  systems. 

4.5.  Given  F  =  d/dt  mVT  ,  F  is  a  force  vector,  m  is  a  constant  mass,  and 
VT  is  the  velocity  vector  of  the  mass  center.  Find  F  ,  F  ,  and  F 

__  _  __  _  x  y  z 

(if  VT  =  Ui  +  Vj  +  Wk  and  w  =  Pi  +  Qj  +  Rk)  with  respect  to  the 
fixed  earth  axis  system. 

4.6.  Given  H  =  Jv  pA  (r  X  V)  dv  where  pAdV  is  the  mass  of  a  particle,  with 

r  as  its  radius  vector  from  the  eg,  and  V  as  its  velocity,  with  respect 
to  the  eg.  Find  Hx  with  respect  to  the  fixed  earth  axis  system. 

4.7.  Write  Hx  in  terms  of  Ix ,  lxy,  and  Ixz,  given: 

Ix  =  Jv  PA  (y2  +  z2 )  dV 

!*y  =  Jv  pa  (*y>  dv 

zxz  =  Jv  Pa  (X2)  dv 


4.8.  Draw  the  three  views  of  a  symmetric  aircraft  and  explain  why  I  =  0, 

■^yz  =  I„  54  0*  What  is  the  aircraft's  plane  of  symmetry? 

4.9.  Using  the  results  from  Problems  4.7  and  4.8,  simplify  the  Hx  equation. 

4.10.  If  G  =  dH/dt  |  XY2 ,  use  the  following: 

H*  -  (plx  -  RIxz  ) 


H  =  QI 

y  y 


H  =  (RI  -  PI  ) 

Z  z  xz  ' 


to  derive  Gx  ,  Gy ,  and  Gz . 


4.11  Define: 

a .  L 

b.  M 

c.  N 

d.  P 


e.  Q 

f.  R 

4.12.  Define  <f>,  0.  what  are  they  used  for?  in  what  sequence  must  they  be 
used?  Explain  the  difference  between  \p  and  8. 

4.13.  What  are  the  expressions  for  p,  Q,  R,  in  terms  of  Euler  angles? 

4.14.  Make  a  chart  that  has  5  columns  and  6  rows.  The  columns  should  contain 
the  terms  of  the  left  hand  side  of  the  equations  of  motion.  Also,  name 
each  equation  (pitch,  drag,  etc.).  Fill  in  all  terms. 

4.15.  What  is  the  difference  between  straight  flight  and  steady  straight 
flight? 
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4.16.  D,  L,  M  =  f(  ,  ,  ,  ,  ) 


Y,  L,  N  ,  ,  ) 


4.17.  Write  L/mU0  in  terms  of  the  stability  parameters.  Define  Lu ,  La,  L., 


L  ,  L.  . 
g '  * 

e 


4.18.  Repeat  4.17  above  for  the  other  5  equations. 


4.19.  Given 


=  m 


[  u  +  Ma+  M»a  +  M..  S  +Mq 

u  a  cc  6  e  q  ^ 


Where 


„  1  3M  „  _  1  3M  „  _  1  3M 

”  "  h  3u  '  “  =  Iy  3«  '  =  Iy  8i 


1_  3M_ 
Iy  3Se 


1  3M 

\  9q 


C  ,  C  ,  C  ,  C  ,  and  C  ; 

jn'm'm#'  m  w  m  9 


(Note:  the  compensating  factor  for  u  is  1/UQ ,  a  and  5e's  compensating 
factor  is  1,  and  a  and  q's  compensating  factor  is  c/2U0 ). 


4.20.  Repeat  Problem  4.19  for 


(Note:  The  compensating  factor  for  3,  Se ,  and  Sr  is  1. 
compensating  factor  for  (3,  p,  and  r  is  b/2U0 ) . 
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4.21.  Given  the  following  representative  wind  tunnel  data: 


De  = 

32.2 

D 

u 

= 

0.00620 

D  = 

a 

0.0000848 

L 

<3 

= 

0 

L 

u 

-  0.0000385 

L 

a 

* 

- 

1.16 

= 

e 

-  0.157 

M 

q 

= 

- 

0.696 

M  = 

u 

0.00104 

M. 

a 

= 

- 

0.000210 

II 

-  18.9 

L. 

a 

= 

0 

M 

a 

-  0.0143 

D6 

= 

0 

e 


Find  a)  A 
e 

b)  G5  (s) 

e 

c)  Gj  (s) 

e 

d)  Na  (s) 

6 

e 
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